Probing dark matter couplings to top and bottom at the LHC 
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Monojet searches are a powerful way to place model-independent constraints on effective oper- 
ators coupling dark matter to the standard model. For operators generated by the exchange of a 
scalar mediator, however, couplings to light quarks are suppressed and the prospect of probing such 
interactions through the inclusive monojet channel at the LHC is limited. We propose dedicated 
searches, focusing on bottom and top quark final states, to constrain this class of operators. We show 
that a search in mono &-jets can significantly improve current limits. The mono-b signal arises partly 
from direct production of 6-quarks in association with dark matter, but the dominant component 
is from top quark pair production in the kinematic regime where one top is boosted. A search for 
tops plus missing energy can strengthen the bounds even more; in this case signal and background 
have very different missing energy distributions. We find an overall improvement by several orders 
of magnitude in the bound on the direct detection cross section for scalar or pseudoscalar couplings. 

PACS numbers: 98.70.Cq 
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I. INTRODUCTION 

Production and detection of dark matter is one of the 
most exciting new physics opportunities at the Large 
Hadron Collider (LHC). The strategy to search for dark 
matter (DM) depends on the physics in the yet-to-be 
fully explored energy range of the LHC. In the maver- 
ick scenario [l[, the DM is the only new particle pro- 
duced and all other new particles arc beyond the scale of 
the LHC. Then the interaction of the DM with standard 
model (SM) particles at these energies can be described 
in terms of an effective field theory (EFT). 

In this case the DM signal at the LHC is missing 
transverse energy (Xr) signals such as mono jets [iHl] or 
monophotons [j, 0] . With an EFT description, one can 
classify all relevant interactions at the LHC in a straight- 
forward way. This scenario also has the advantage that 
the connection between DM annihilation, direct detec- 
tion, and collider signals is simple. 

The ATLAS [1, Q and CMS [H [HI collaborations 
have published monojet constraints on the scale of new 
interactions in the EFT, which are then used to place 
constraints on the DM-nucleon scattering cross section. 
These constraints are most effective for DM masses be- 
low 100 GeV. Meanwhile, there has been rapid progress 
in the direct detection of dark matter [Tj, II^Ij with 
the strongest bounds on DM-nucleon scattering at DM 
mass of around 50 GeV and for spin-independent scat- 
tering. These two approaches are complementary, and 
connec ting them has been the focus of many recent stud- 
ies [11, [If. 



challenging to constrain the scalar operator, where inter- 
actions between dark matter and quarks are mediated by 
a heavy scalar mediator: 
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summing over all quarks.^ The form of the interac- 
tion is fixed by minimal flavor violation (MFV) [T6| . 
Scalar interactions with SM quarks are typically strongly 
constrained by flavor changing neutral current measure- 
ments, but in MFV these dangerous flavor violating ef- 
fects are automatically suppressed. 

Because the interactions arc proportional to quark 
mass, however, the monojet4>E^T signal rate appears to 
be suppressed by the light quark masses. The ATLAS 
monojet search based on 4.7 fb^^ at 7 TeV sets a limit 
of > 30 GeV [1], including only couplings to charm 
and lighter quarks. This bound is much weaker than 
constraints on operators mediated by vector or axial in- 
teractions. 

In this paper, we point out that the direct search for 
production of dark matter in association with third gen- 
eration quarks can enhance the reach of the LHC for 
dark matter coupled to quarks through a scalar interac- 
tion. Direct b production gives rise to a mono fe-jet sig- 
nal. We also show that the kinematics of top quark pairs 
plus dark matter is such that boosted tops may form the 
dominant contribution to the mono-6 signal. However, 
monojet searches veto on more than 2 hard jets, so a 



A. Scalar Operator 

While the monojet search is extremely effective for 
many of the possible operators, it is not necessarily the 
optimal way to study all of them. In particular, it is 



There are closely related operators, for example the pseudoscalar 
operator ■j^q-y^qX'y^ X , which will have almost identical collider 

constraints. The direct detection cross section for these operators 
is velocity suppressed, however, so the best limits will come from 
the LHC. It should be understood that our limits apply to both 
scalar and pseudoscalar operators. 



FIG. 1: Some of the dominant diagrams contributing to associated production of DM with bottom and top quarks. 



better strategy to probe the couplings to top quarks is 
the study of ti + final states. 

The scalar interaction has also been studied recently 
in Ref. [TtI , which showed that heavy quarks in loops can 
significantly enhance inclusive monojet production. We 
focus instead on direct identification of the heavy quarks 
in the final state. Mono-6 final states from dark matter 
have also been studied in Refs. [H, [3, although not in 
the context of MFV, so the top quark contribution to the 
mono-5 signal was not considered. 

In Section [IT] we study the mono fe-jct signal where the 
leading jet is 6-tagged. This search can improve con- 
straints on the DM-nucleon cross section, cr„, by several 
orders of magnitude compared to current ATLAS limits. 
In Section [IIII we show an even stronger limit can be ob- 
tained from a search for dark matter in association with 
top quarks, tt+^r- This is also the final state studied in 
searches for stops, supersymmetric partners to tops, and 
we use published results to derive limits. We find that 
the limit on cr„ is stronger by another factor of approxi- 
mately 2 compared to the mono &-jet search. 



II. MONO &-JET SEARCH 

The scalar operator gives rise to 6-jets plus 0t via 
direct b production, as well as from production of top 
quarks which then decay. Direct h production occurs 
through h and gluon-initiated processes, such as bg — > 
XX -\- b; several example diagrams are shown in Fig. [T] 
In comparison to the light quark initial states, these 
processes are suppressed by the 5-quark parton density. 
However, the enhancement due to the MFV form for the 
coupling is more than enough to compensate this. 

Furthermore, gg — )■ XX + ti turns out to be the dom- 
inant contribution to the monojet signal. Thus, the fi- 
nal states are highly 6-enriched. At the same time, fo- 
cusing on exclusive 6-taggcd final state reduces the SM 
backgrounds significantly. Therefore, we expect an im- 
provement in the LHC reach for the scalar operator by 
requiring a 6-tagged monojet. 

Before presenting our results, we summarize our event 
simulation methods. We use MadGraph 5 [20| for parton- 
level cross sections, interfaced to Pythia 6 [2l| for show- 



ering and hadronization, and Delphes 2 [2^ for detector 
simulation. For Delphes, we set a 60% tagging efficiency 
for 6, 10% mistag for c, and 0.2% mistag rate for light 
quarks and gluons [1^. Jets are clustered into R ~ 0.4 
anti-fc^ jets. 

Up to two hard jets are allowed in the monojet and 
mono-6 searches, so we must consider Next-to-Leading- 
Order (NLO) corrections in our simulation of the signal. 
We generate matched samples with fcr-jet MLM match- 
ing. For SM backgrounds, we generate W/Z and ti with 
up to 2 jets. For the signal, we generate XX+iets, in- 
cluding up to 2 jets, for all fiavors other than tops. We 
separately include XX + ti at leading order. Finally, we 
normalize all matched samples with NLO cross sections 
computed using MCFM [Ij] . 

For the signal region we require 0t > 350 GcV, a 
leading 6-taggcd jet with pT > 100 GcV, |?7| < 2.5, and 
no isolated leptons. We also allow an additional softer 
jet. but no more than two jets with px > 50 GcV. There 
is a cut on the azimuthal separation between and 
the second jet, A(j){^T,PT) > 0-4, in order to suppress 
the mismeasured dijet background. This signal region 
overlaps well with those used in previous studies, and 
furthermore the dependence on the cut values appears to 
be mild. 

The resulting cross sections at 8 TeV are given in Ta- 
ble m We have split the signal into three contributions: 
coupling to charm and light quarks, direct b and bb pro- 
duction from coupling to b, and ti production. 

Associated production of DM with ti constitutes the 
dominant signal for both monojet and mono-5 signals 
because of the enhancement from the top mass and be- 
cause of the production of boosted tops which can be 
tagged as 6-jets. Events where only one top is boosted, 
and where the other top gives rise to low pT jets, can 
pass the niono-6 cuts. Of the events that pass the 
requirement and lepton veto, 17% of events survive the 
veto based on the px of the third jet. In comparison, 
80% of events from direct 6 production survive the jet 
veto. In both b and top production, about 50% of those 
events then have a leading jet which is 6-tagged. Note 
that this assumes the same 6-tagging efficiencies for the 
6-jets inside the boosted tops. Without dedicated study 
by the experimental collaborations, this is an idealized 
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Process 


Mono jet 


6-tag 


6-tag on ji 




Z+jets(fake) 


406 fb 


11 fb 


7 fb 




Z+b+jet 


6.7 fb 


4 fb 


3 fb 


Background 










W+jets,W^+b 


95 fb 


3 fb 


2 fb 




tt+jets 


16 fb 


11 fb 


6 fb 




XX+jeta 


11 fb 


0.9 fb 


0.7 fb 


Signal 


XX + 6+jets 


65 fb 


40 fb 


33 fb 




XX + tt 


244 fb 


156 fb 


113 fb 



8 TcV, 20 lb-\ PT(b-jct) > 100 GcV 



TABLE I: Monojet and mono-6 search at 8 TeV: Cross 
sections for dominant backgrounds and signal with cuts of 
Xt > 350 GeV, > 100 GeV as described in the text. 
For the signal we take M* = 50 GeV and mx = 10 GeV. 
The row XX+ jets includes only DM coupling to charm and 
lighter quarks. Note that XX + b+ jets includes single &-jet 
and bb production. In the column labeled 6-tag, a 6-tag on 
any jet withpT > 50 GeV is required, while in the last column 
the leading jet must be 6-tagged; this choice does not lead to 
significantly different results in setting limits. 



assumption. Wc note that most of the top jets are from 
mildly boosted tops, with have pT ~ 400 GcV. Hence, 
wc do not expect the 6-tagging efficiency to degrade sig- 
nificantly. At the same time, it may be possible to use 
additional information to tag these jets as coming from 
boosted tops. We will leave this for a future study. 

The dominant SM backgrounds are Z{i'D)+jets and 
W-|-jets, where a jet fakes a &-jet. Although these are 
suppressed by the mistag rate, they are still a larger 
contribution than direct b production from Z + b and 
tt backgrounds. Depending on the 6-tagging algorithms 
used, however, it may be possible to further reduce the 
background from Z/W+jets. 

Kinematic distributions in JS^t and leading jet pT are 
shown in Fig. [2l The J^t spectrum for b production is 
very similar to that for Z-|-jcts. despite the fact that 
the signal arises from a contact interaction. This is 
partly because in the signal case, the initial states in- 
clude sea quarks, while for the dominant Z^-l-jets back- 
ground roughly 70% of events are initiated by at least 
one valence quark. Meanwhile, XX + tt final states tend 
to have more J^t because of the requirement of producing 
massive top quarks. 

To estimate the expected bound, we compute the num- 
ber of signal events such that < 2.71 to obtain a 
90% CL bound [1]. A systematic uncertainty of 5% 
is assumed. We also compute bounds for 14 TeV and 
100 fb~^, keeping the same cuts. A higher cut can 
improve bounds in the inclusive monojet case @ ; in this 
case, however, stronger cuts would also require higher px 
6-jets, where the 6-tagging efficiency can degrade. 

Fig.[3]shows our constraints on M^, in the left panel and 
the corresponding limits for direct detection in the right 
panel. The scalar operator gives rise to spin-independent 
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FIG. 2: Mono-6 search at 8 TeV: Distributions for^T and 
, the transverse momentum of the leading 6-jet, for some 
of the dominant SM backgrounds and the signal. We separate 
the DM signal into contributions from direct 6 production and 
from ti production. For the signal we take M, = 50 GeV and 
mx = 10 GeV. 



DM-nucleon scattering, with a cross section of 



^ (Q-^^- j^^^ . 2 X 10-3W r^V, (2) 



7tM2 



A/* J 



compared to XENONIOO limits of about IQ-'^^cm^ at 
mx = 10 GeV We find a factor of 8 improvement 
in cr„ limits with a mono-& search compared to an in- 
clusive monojet search. Overall, we find a factor of 250 
improvement compared to an inclusive search where only 
the coupling to charm quarks and lighter is considered (as 
in the most recent ATLAS study Q). 
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FIG. 3: (Left) Expected 90% CL limits on the scalar operator from a mono-b search, including couplings to tops and bottoms. 
For the mono-6 search at 8 TeV we also show the limit if b-jets from top production are not included (dotted line). For 
comparison we include limits for an inclusive monojet search with no 6-tag, and current ATLAS limits from §]. (Right) 
Corresponding constraints on the spin- independent nucleon scattering cross section, along with XENON 100 limits [12], and 
projected sensitivity for XENONIT 



III. TOPS PLUS MISSING ENERGY SEARCH 

As shown in the previous section, the process gg 
XX + tt contributes the dominant component of the 
monojet and the mono-6 signals. The monojet and mono- 
h searches veto on more than two high-pT jets, however, 
cutting out a large fraction of ti events. A stronger con- 
straint on this couphng can be obtained from dedicated 
searches. 

Models of supersymmetry also have a signature of top 
pairs plus missing transverse energy. We apply the recent 
ATLAS 8 TcV search for top-quark superpartners with 
1-lepton final states [1^ using 13 fb~^ of data to these 
scalar dark-matter couplings.^ The signal regions require 
1 isolated lepton,^T > 150 GeV, transverse mass^ my > 
120 GeV, 4 jets with pr > (80,60,40,25) GeV and at 
least 1 6-tag. 

Fig. m shows the and rriT distributions of the sig- 
nal and the dominant background, tt. The DM signal 
is significantly harder in the spectrum, whereas the 
background is highly peaked towards low J^t because the 
primary source of is from the neutrinos in the top de- 
cay. Meanwhile, it is unlikely that stronger cuts on niT 
above 120 GeV would substantially improve the ratio of 
signal to background. 



^ We have also calculated constraints using the CMS 1-lepton fi- 
nal state search [27l | and obtain limits that are similar although 
slightly weaker. 

^ The transverse mass is defined as (my)-^ = 2p'^^^q^{l — cos A</<) 
with A(p the azimuthal separation between lepton and missing 
momentum directions. 



We find the best constraints come from the signal re- 
gion D (SRD) of the ATLAS study, wither > 225 GeV. 
Although there is another signal region with J^fr > 275 
GeV, the systematic uncertainties increase significantly. 
We thus apply the ATLAS SRD cuts to simulated data 
to derive our 13 fb""'^ limits. The signal cross section with 
these cuts is 

O-signal = 173 fb (3) 

assuming M^, = 50 GeV and mx = 10 GeV. We give 
limits on and (t„ in Fig. O Because uncertainties are 
systematics dominated for this signal region, we do not 
expect a significant improvement of limits with 20 fb^^ 
of data. 

Fig. El also shows limits for 14 TeV with 100 fb"^ 
of data, keeping the same cuts as above. We simulate 
VF-|-jets in addition to tt for our background estimate, 
and assume the systematic error on the background is 
the same as in the 8 TeV analysis. We also calculated 
constraints for a search with an all-hadronic final state 
[28l [29| ; in this case it may be possible to improve the 
bounds on M, by 10-20%, depending on the detector ac- 
ceptances and systematic uncertainties. 

IV. DISCUSSION 

We have shown that limits on scalar (and pseudoscalar) 
interactions of dark matter with quarks can be improved 
significantly by directly searching for final states with b- 
jets and tops. Compared to an analysis including only 
light quarks, we find a factor of 400 improvement in lim- 
its on (T„, and compared to an inclusive monojet search 
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FIG. 4: {LeftjJS'T distribution after requiring an isolated lepton and rriT > 120 GeV. (Right) Transverse mass rriT distribution 
requiring an isolated lepton and J^t > 150 GeV. The dark matter mass is mx ~ 10 GeV. 
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FIG. 5: (Left) Expected 90% CL limits on scalar operator from applying a search for supersymmetric tops with one lepton in the 
final state. 8 TeV limits are obtained using the results of [1^. Also shown are ATLAS limits from 8!|. (Right) Corresponding 
constraints on nucleon scattering cross section, along with XENONIOO limits [l2|, and projected sensitivity for XENONIT i25| . 



including couplings to all quarks, we find a factor of 15 
improvement. For 8 TeV data, the corresponding con- 
straints on direct detection are below the regions favored 
for light dark matter interpretations of DAIVIA [s^l and 
CoGeNT [m. 

Couplings to heavy quarks can also lead to an enhance- 
ment of inclusive monojet production through loops; 
Ref. O found > U8t\l GeV for small DIVI mass 
using 7 TeV data. However, these loop corrections as- 
sume that the operator is generated by a heavy neutral 
scalar. Although our constraints are weaker, the searches 
discussed here directly probe couplings of dark matter to 
top and bottom. Furthermore, the J^^t spectrum for the 
tt final state is strikingly different from the background. 



It may be possible to use the difference in shapes to im- 
prove limits from the searches discussed here. 



Finally, in this paper we have assumed a contact inter- 
action for simplicity. As discussed in Refs. this 
assumption must be compared to the derived bounds on 
A'/*. In this case, the best limit we obtain at 8 TeV is 
A'/* > 110 GeV, and for this value a significant fraction 
of events (over 50%) violate the criteria in Refs. jsl. IT3|. 
A UV completion for this operator is necessary to derive 
fully consistent constraints. At the same time, the re- 
sults will be more model-dependent and we reserve this 
analysis for future work. 
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